Localization of fermions in an anisotropic random magnetic field 
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We study the localization of fermions in an anisotropic random magnetic field in two dimensions. It 
is assumed that the randomness in a particular direction is stronger than those in the other directions. 
We consider a network model of zero field contours, where there are two types of randomness - the 
random tunneling matrix element at the saddle points and unidirectional random variation of the 
number of fermionic states following zero field contours. After averaging over the random complex 
tunneling amplitude, the problem is mapped to an SU(2N) random exchange quantum spin chain 
in the N — > limit. We suggest that the fermionic state becomes critical in an anisotropic fashion. 
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The localization problem of a nonrelativistic particle in 
a random perpendicular magnetic field has been a sub- 
ject of intensive research and debate. In particular, the 
problem arises from the study of disordered half-filled 
Landau level in the quantum Hall regime §. In the ab- 
sence of disorder, the Chern-Simons theory of electrons in 
the half-filled Landau level leads to an effective theory of 
composite fermions in the fluctuating gauge field |jj . Due 
to the fact that the density distribution of the fermions 
are tied to the fluctuating magnetic field, static disorder 
leads naturally to the random static magnetic field with 
zero mean Since the longitudinal resistivity of the 
electrons in the half-filled Landau level is the same as 
that of composite fermions, the problem of localization 
of composite fermions in a random static magnetic field 
is important for understanding transport properties at 
the half-filling. 

One might expect that the system belongs to the uni- 
tary universality class as the disordered system with bro- 
ken time reversal symmetry and zero Hall conductance. 
According to conventional scaling theory, all the states 
are localized in two dimensions for the unitary class ||. 
It is established that those systems can be described by a 
unitary non-linear sigma model (NLcrM) without topolog- 
ical term in the long wavelength limit. Some years ago, 
Zhang and Arovas j| suggested that local fluctuation of 
induce the long-range interactions between topo- 
logical densities in the NLcrM. It was claimed that this 
may give rise to a novel derealization transition when 
the conductance at short distances becomes larger than 
a critical value of the order of e 2 /h. 

Since then, it has been controversial whether the model 
has a delocalized state or not. There have been several 
numerical and analytical studies leading to different con- 
clusions. Some of them conclude that the system indeed 
belongs to the unitary class and all the states are local- 
ized ^pQ. On the other hand, some numerical calcula- 
tions on the lattice model reached the conclusion that 
there are delocalized states near the band center [0,0. 



The very recent extensive numerical calculation by Fu- 
rusaki |)| showed that a state in the band center is not 
localized when there is a special particle-hole symmetry 
in the lattice model. On the other hand, in the absence 
of the particle-hole symmetry, all the states are local- 
ized. This suggests that there is no extended state in the 
generic continuum model j9|. 

Previously Furusaki, Lee, and the author |6| map a 
network model || of snake states (which are fermionic 
states following zero field contours) to a quantum spin 
chain. They suggest that the quantum spin chain has a 
gap in the excitation spectrum so that the fermions in 
the original problem is always localized. In this paper, 
we elaborate this mapping further and point out that 
some ideas of Zhang- Arovas theory may be realized in an 
anisotropic version of the random flux problem. In par- 
ticular, we suggest that the number of snake states can 
fluctuate in space, which was set to be a fixed number in 
previous numerical and analytical studies. We show that 
if the number of snake states changes only in a partic- 
ular direction, the fermionic state becomes critical and 
the critical behavior could be different for different direc- 
tions. 

Let us first review the mapping from a network model 
of snake states to a quantum spin chain pjl0||. We as- 
sume that the random magnetic field has zero mean and a 
non-zero variance. We consider the semiclassical limit in 
which the correlation length of the random magnetic field 
is large compared to the typical magnetic length. In this 
limit of smooth disorder, a fermion moves along contours 
of constant magnetic field. The contours which percolate 
geometrically across the system would be B = contours. 
It turns out that semiclassical motion of a fermion along 
B = contours follows snake-like trajectories. In the net- 
work model, these zero field contours are represented as 
links on a square network, carrying the fermionic "snake" 
states The node of the network corresponds to the 
saddle points of the magnetic field where the scattering 
to either the right or left occurs B. The quantum inter- 
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ference due to these scattering processes are modeled by 
random tunneling amplitudes at the saddle points and 
the mixing of the snake states. 

It is important to notice that the fermionic states along 
the zero field contours arise in pairs at a given energy and 
they are propagating in the same direction. This is be- 
cause these states are generated from a symmetric and 
antisymmetric linear combinations of the Landau levels 
on either side of the contour. This implies that there 
exist two snake states per each local Landau level be- 
low the Fermi energy. Therefore, there are always even 
number, M, of snake states in each link of the network 
model. Previous studies assumed that the number of 
snake states, M, is fixed HI). However, notice that M 
could be a spatially varying quantity as far as this is an 
even number because the magnitude of the magnetic field 
changes in space. Later, we argue that if the change of M 
occurs only in one direction, then it leads to a significant 
effect. At the moment, we assume that M is a fixed even 
number. 

In the network model, two types of random com- 
plex tunneling amplitude are introduced, tj and tjk- tj 
represents the tunneling between the snake states, ipj 
(j = 1, M), in neighboring zero field contours and tjk 
corresponds to the mixing between snake states in differ- 
ent channels (among M number of channels) in the same 
zero field contour. Using the replica trick, the disorder 
averages over tj and tjk are done. Taking y as the imagi- 
nary time, one can regard the disorder averaged action as 
representing interacting fermions, tp a j (j = X,...,M and 
a = 1, 27V) on M number of coupled one dimensional 
systems. Here TV is the replica index. Substituting the 
fermion operators with the generators of SU(2N): 

S tf3 = V4" i>/3j " <W E ^7 j > ( 1 ) 
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we can obtain M coupled SU(2N) spin chains given by 
the Hamiltonian ||: 

ff = £X>Tr[^(z + 

X j 

X j<k 

where a is the lattice spacing corresponding to typi- 
cal distance between snake snakes. Here Jj — (\tj\ 2 ), 
Jjk = -(\tjk\ 2 ), and Tr[AB] = Y. a p A l B % any 
2N x 2N matrix A and B. In the case of decoupled 
chains, each spin chain is described by a particular rep- 
resentation of SU (2 AT), characterized by a Young tableau 
with a single column of length N . 

Notice that there are antiferromagnetic couplings 
(Jj > 0) in each spin chain j and ferromagnetic cou- 
plings {Jjk < 0) between the spin chain j and k for 
each x. The spin at each site represents a snake state 



at a given time at that particular position. The antifer- 
romagnetic coupling, Jj, comes from the fact that the 
tunneling between the snake states in the same chan- 
nel always occurs between oppositely propagating snake 
states. On the other hand, the mixing between snake 
states in different channels happens for snake states prop- 
agating in the same direction, which leads to the ferro- 
magnetic coupling, Jjk- Notice that each chain corre- 
sponds to a single-channel network model with non-zero 
Hall conductance. In order to impose the constraint of 
zero Hall conductance we require the alternation of the 
bond strength on the neighboring chains to be staggered. 
This means that Jj(x = even) = Jj + \{x = odd) and 
Jj(x — odd) = Jj + i(x = even). In terms of the original 
problem, the localization of the fermions corresponds to 
the existence of an energy gap in the spin chain (in the 
replica limit, N — > 0). 

Due to the fact that there are M number of snake 
states in each link, the conductance of the system at 
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short distances is G — Mv, Thus the large conduc- 
tance limit, M 3> 2, corresponds to the large number of 
coupled chains. The interchain coupling is ferromagnetic 
and it exists for any two of M spin chains. Thus, as M 
becomes large we expect that, the coupled chains become 
essentially single spin chain with a uniform antiferromag- 
netic coupling J 0: 

H = J2 J^[S(x + a)S(x)} , (3) 

X 

where 5 (x) is in the totally symmetric representation de- 
scribed by a Young tableau with AI number of columns 
of length N. It is useful to think about the SU{2) case 
(N =1). In this case, M corresponds to 25, where 5 is 
the ordinary spin quantum number. Since M is always 
an even number, the system can be described by an inte- 
ger spin chain. Therefore, there is a gap in the excitation 
spectrum in the spin chain. If one assumes that N — > 
limit is not singular, this may imply that the electrons 
are localized in the original problem even for large bare 
conductances. 

Now let us consider the case that the number of snake 
states is not fixed, but can change in space. In particular, 
we assume that the number of snake states changes only 
in one direction. We will discuss the isotropic case later. 
In this anisotropic random flux problem, we can take the 
x direction as the direction in which M changes in space. 
We assume that there are two types of randomness in the 
problem - one of them is the random tunneling amplitude 
at the saddle points and the other one is the number of 
snake states. After performing the average over the ran- 
dom tunneling amplitude, we get 

ff = E E E J]Tr[Sj(x + a)Sj(x)} 

V l p <x<r p l<j<M p 

+E E jjnSj(i P+ i)Sj(r p )] 

P l<j<Min(Mp,M p+ i) 
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+ E E E JjkTr[Sj(x)S k (x)} . (4) 

p l p <x<r p l<j<k<M p 

Here the coupled spin chains consist of many subsystems 
labeled by p. l p and r p are the positions of the left and 
right ends of the pth subsystem (l p +i = r p + a). Each 
subsystem labeled by p has M p number of coupled spin 
chains. In general, M p+ \ — M p = M p is an even number. 
In the large bare conductance limit, where M p is much 
larger than 2, the effective Hamiltonian can be written 
as 

ff = E E J P ^[Sp(x + a)S p (x)] 

+ Y,Mm(J p ,J p+1 )Tv[S p+1 (l p+1 )S p (r p )} . (5) 
p 

Here the system still consists of many subsystems. Each 
subsystem labeled by p is a uniform antiferromagnetic 
spin chain in a representation described by a Young 
tableau with M p number of columns of length N. 

It is instructive to think about SU{2) case (N = 2). In 
this case, the above Hamiltonian represents the coupled 
antiferromagnetic integer spin chains. Each integer spin 
chain has a uniform antiferromagnetic coupling J p and 
the spin magnitude S p = M p /2. Thus each subsystem 
would have a gap in the excitation spectrum if it were 
an infinite system. Notice that the second term in the 
Hamiltonian corresponds to the coupling between edge 
spins of adjacent spin chains. For an integer spin chain 
with finite size and spin magnitude S Pl the edge spin has 
the spin magnitude S p /2 fPlJ| , This can be seen from the 
Berry phase contribution to the effective action for an 
integer spin chain of length L: 

S BP = S p J2(-l) x/a n(x) ~ f / dx^ 

X J 

= ^[[O(L)-n(0)]+47rQ] , (6) 

where Vl(x) is the solid angle subtended by the closed 
path on the surface of a unit sphere defined by the time 
evolution of a unit vector at x. Q is an integer measur- 
ing the number of times the spin configuration covers the 
surface of the sphere. Since % x inQ — 2ttQ, the last 
term has no effect and can be dropped. 

At the interface x = L, between the spin chains with 
S p and S p+ i, the Berry phase contribution in the limit 
L ^> a becomes 

s BP = ^n(L) - ^n(L + a) « Sp ~ 2 Sp+1 n(L) . 

(7) 

Since (S p — S p+ \)/2 = (M p — M p+ i)/A, the edge spin 
can be both the half-integer and integer spin. In the sys- 
tem described by Eq.(pf), this generates edge spins with 



random spin magnitude Sp dse (X) at the interface, X, 
between adjacent pth and (p + l)th integer spin chains. 
Furthermore, depending on whether r p — l p is even or 
odd, the residual couplings between those edge spins can 
be either antiferromagnetic or ferromagnetic. The mag- 
nitude of the coupling goes as J~ p e~( rp ~ lp " l ^ p , where £ p 
is the correlation length. From the above consideration, 
one can see that the effective coarse-grained Hamiltonian 
of the system in the case of SU(2) is given by 

H = J2 J{X)S cdsc (X + A) ■ S cdgc (Y) , (8) 
x 

where J has a distribution of both antiferromagnetic 
{J > 0) and ferromagnetic couplings (J < 0). Here the 
spin magnitude can be also random in space and A is the 
typical distance between interfaces of different subchains. 

Coming back to SU (2N) case, one can generalize the 
above argument. For example, at the interface, X, be- 
tween two adjacent spin chains, the edge spin is gener- 
ated and it is in the representation characterized by a 
Young tableau with M cdg0 = (even number) /2 number 
of columns of length N. At the same time, the sign and 
magnitude of the residual spin exchange interaction be- 
tween those edge spins can be random. Thus the effective 
Hamiltonian can be written as 

H = J2 J{X)Tr[S cdsc (X + A)S cdsc (X)} . (9) 
x 

The number of columns of the Young tableau for the edge 
spins at X has a random distribution. 

It is rather easy to understand what the effective 
Hamiltonian stands for in the original snake-states pic- 
ture. In the region where the number of snake states 
changes from one even number to the other, the num- 
ber of remaining snake states after forming all possible 
localized configuration out of available snake states, is 
still always even and they are all in the same direction. 
Since most of the snake states in the bulk are localized, 
the remaining degrees of freedom are described by the 
weak coupling across a long distance between residual 
snake states at the interface between two region of dif- 
ferent number of snake states. The number of residual 
snake states can be an arbitrary even number. In the 
limit of strong mixing, the residual snake states can be 
represented by a single large snake state. In the spin 
chain language, the residual snake states correspond to 
S' edge . Thus the number of columns in the Young tableau 
corresponds to the number of residual snake states at the 
interface. If the weak coupling occurs between the resid- 
ual snake states propagating in the same direction, the 
coupling in the spin chain is ferromagnetic. For the snake 
states propagating in the opposite direction, the coupling 
becomes antiferromagnetic. 

The SU(2) case described by the Hamiltonian Eq.(|^) 
was studied by various methods 0,0]. It turns out 
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that the real space renormalization group (RG) [ [l4||T5| 
treatment can provide a lot of insight about the ground 
state of the system |L2l . Let us consider an isolated link 
with two sites, where two neighboring spins are coupled 
by J. The ground state of an isolated two-site system 
would have maximum (jj < 0) or minimum (jj > 0) spin 
quantum number depending on the sign of the exchange 
coupling. There would be an energy gap A to the first 
excited state. Since the total system consists of these 
links, the system can be described by a distribution of 
energy gaps and the spin magnitude. In the real space 
RG scheme, one successively replaces the strongest links, 
where the energy gap is the largest, by effective spins. 
The RG transformation preserves the form of the Hamil- 
tonian, but it changes the probability distribution of the 
gaps (or exchange couplings) and the spin magnitude. 

In the low energy limit, the RG flows reach a single uni- 
versal fixed point and the distributions become universal 
at the fixed point [Q. At low energies, the random spin 
chain can be described by weakly coupled large effective 
spins [ p^[ . As the links are replaced by effective spins, 
the effective couplings or the gaps of the links decrease. 
At the same time, the average distance between neigh- 
boring spins and the magnitude of the effective spins in- 
crease. It was found that the average size of the effective 
spins diverges as T~ a while the average distance between 
neighboring effective spins increases as T~ 2a [Q. That 
is, the correlation length diverges in the zero tempera- 
ture limit. This suggests that the following correlation 
function decays as a power law at the zero temperature 

& 

c(x -y) = (£ XY s{x) • s(y)> oc lx * Yl „ , (io) 

where £ XY = n^sgn(- J{Z)). 

In the original problem of fcrmions in an anisotropic 
random flux distribution, this may imply that the 
fcrmionic state is critical in the direction of changing 
number of snake states or large variation of the mag- 
nitude of random magnetic field. This means that the 
wavefunction is not exponentially localized in this direc- 
tion, but decays as a power law as a function of system 
size. Therefore, in some sense, it behaves like a critical 
state at the transition between metal and insulator. The 
reason why the correlation becomes critical is that some 
rare events across the system in the x-direction domi- 
nates the disorder averaged correlation function. Since 
the spin chain does not have an excitation gap, the dis- 
order averaged correlation in the time direction cannot 
decay exponentially. Due to the critical nature of the 
spin chain, we expect that it decays also as a power law 
even though the power law would be in general different 
for different directions. Thus, the fermionic state in the 
y-direction is also critical and the critical behavior may 
be different from that in the x-dircction. In the absence 



of detailed information about the exponents v for differ- 
ent directions, it is hard to tell in which direction the 
wavefunction decays faster. 

In the isotropic problem, the number of snake states 
can change in any direction. In the spin chain language, 
it will induce the time and space dependent random spin 
size and the couplings. Thus one gets an annealed disor- 
der problem instead of a quenched one. Even though it is 
not clear what the effects of this type of disorder would 
be, it is possible that the effects of annealed disorder 
could be much weaker so that the random-tunncling-only 
model may be enough to describe the isotropic random 
flux problem. 

In conclusion, we consider the motion of fermions in an 
anisotropic random magnetic field. We assume that the 
number of fcrmionic states following zero field contours 
can be changed only in one direction. The mapping from 
a corresponding network model to a random quantum 
spin leads to the conclusion that the fcrmionic states are 
critical in an anisotropic fashion. 
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